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Chern-Simons (CS) modified gravity is an extension to general relativity (GR) in which the metric 
is coupled to a scalar field, resulting in modified Einstein field equations. In the dynamical theory, 
the scalar field is itself sourced by the Pontryagin density of the space-time. In this paper, the 
coupled system of equations for the metric and the scalar field is solved numerically for slowly- 
rotating neutron stars described with realistic equations of state and for slowly-rotating black holes. 
An analytic solution for a constant-density nonrelativistic object is also presented. It is shown that 
the black hole solution cannot be used to describe the exterior spacetime of a star as was previously 
assumed. In addition, whereas previous analysis were limited to the small-coupling regime, this 
paper considers arbitrarily large coupling strengths. It is found that the CS modification leads to 
two effects on the gravitomagnetic sector of the metric: (i) Near the surface of a star or the horizon 
of a black hole, the magnitude of the gravitomagnetic potential is decreased and frame-dragging 
effects are reduced in comparison to GR. {ii) In the case of a star, the angular momentum J, as 
measured by distant observers, is enhanced in CS gravity as compared to standard GR. For a large 
coupling strength, the near-zone frame-dragging effects become significantly screened, whereas the 
far-zone enhancement saturate at a maximum value AJmax ~ {AI / R)Jgr- Using measurements of 
frame-dragging effects around the Earth by Gravity Probe B and the LAGEOS satellites, a weak 
but robust constraint is set to the characteristic CS lengthscale, < 10* km. 



I. INTRODUCTION 

Gravity, one of the four fundamental forces of nature, 
is elegantly described by Einstein's theory of general rel- 
ativity (GR). Nearly a century after its discovery, GR 
has successfully passed the more and more subtle and 
precise tests that it has been submitted to (for a review, 
see for example Ref. [T]). Nevertheless, Einstein's theory 
is probably not the final word on gravity. We expect that 
a more fundamental theory unifying all forces should be 
able to describe not only gravity but also quantum phe- 
nomena that may take place, for example, at the center 
of black holes. Since gravity does seem to describe nature 
quite faithfully at the energy and length scales that are 
accessible to us, we expect that it may be the low-energy 
limit of such a fundamental theory. If this is the case, 
gravity should be described by an effective theory, for 
which the action contains higher-order curvature terms 
than standard GR, the effect of which can become ap- 
parent in strong gravity situations. 

One such theory is modified Chern-Simons (CS) grav- 
ity (for a review on the subject, see Ref. |2]). In this 
theory, the metric is coupled to a scalar field through the 
Pontryagin density RR (to be defined below). The dy- 
namical and nondynamical versions of the theory (which 
in fact are two separate classes of theories), then differ in 
the prescription for the scalar field. In nondynamical CS 
gravity, the scalar field is assumed to be externally pre- 
scribed. It is often taken to be a linear function of coordi- 
nate time (the so-called "canonical choice" ) , which selects 
a particular direction for the flow of time |3j , and induces 
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parity violation in the theory. Nondynamical CS grav- 
ity then depends on a single free parameter, which has 
been constrained with measurement of frame-dragging 
on bodies orbiting the Earth 4 , and with the double- 
binary-pulsar O E]. Nondynamical CS theory is quite 
contrived as a valid solution for the spacetime must sat- 
isfy the Pontryagin constraint RR = 0; it should there- 
fore rather be taken as a toy model used to gain some 
insight in parity-violating gravitational theories. 

Dynamical CS gravity, which is the subject of the 
present paper, is a more natural theory where the scalar 
field itself is given dynamics (even though arbitrariness 
remains in the choice of the potential for the scalar field) . 
The scalar field evolution equation is sourced by the Pon- 
tryagin density, which is non-vanishing only for space- 
times which are not refiection-invariant, as is the case 
in the vicinity of rotating bodies. Contrary to the non- 
dynamical theory, though, dynamical CS gravity is not 
parity breaking, but simply has different solutions than 
GR for spacetimes which are not refiection-invariant (see 
discussion in Sec. 2.4 of Ref. [2 ). Dynamical CS gravity 
has only recently received some attention, as it is more 
complex than the nondynamical version. Refs. [TjlSl com- 
puted the CS correction to the Kerr metric in the slow- 
rotation approximation. Ref. [S| studied the effect on CS 
gravity on the waveforms of extreme- and intermediate- 
mass ratio inspirals. Recently, Ref. [101, proposed a so- 
lution for the spacetime inside slowly-rotating neutron 
stars, assuming that the solution outside the star was 
identical to that of a black hole of the same mass and an- 
gular momentum. All the aforementioned studies were 
done in the small-coupling limit, i.e. considering the CS 
modification as a perturbation around standard GR. 

The two main points of the present work are as fol- 
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lows. First, we show that, in contrast with GR, the space- 
time around a slowly-rotating relativistic star is different 
from that of a slowly-spinning black hole in CS gravity, 
owing to different boundary conditions (regularity con- 
ditions at the horizon for a black hole versus continu- 
ity and smoothness conditions at the surface of a star). 
Second, we solve for the CS modification to the metric 
and the scalar field simultaneously, in the fully-coupled 
case (nonperturbative with respect to the CS coupling 
strength), for a slowly-rotating star or black hole. Our 
motivation in doing so is that frame-dragging effects are 
difficult to measure and are not highly constrained; it is 
therefore still possible that they differ significantly from 
the GR prescription. The solutions obtained would be 
exact (modulo the slow-rotation approximation) if CS 
gravity were taken as an exact theory. If one asumes the 
CS action is only the truncated series expansion of an 
exact theory, then our solutions are only meaningful at 
the linear order in the CS coupling strength. 

As shown in previous works, the CS correction only 
affects the gravitomagnetic sector of the metric at lead- 
ing order in the slow-rotation limit. We find that for 
both stars and black holes, the CS correction leads to a 
suppression of frame-dragging at a distance of a few stel- 
lar radii, or a few times the black hole horizon radius. 
This suppression is perturbative in the small-coupling 
regime but can become arbitrarily large in the nonlinear 
regime, asymptotically leading to a complete screening 
of frame-dragging effects near the star or black hole in 
the large-coupling-strength limit. Far from the bound- 
ary of a star, the magnitude of frame-dragging effects is 
enhanced: we find a correction to the t(j) metric compo- 
nent Ag*'^ sa -2AJcs/r^ at large radn, with AJqs > 0. 
This means that for a given angular rotation rate 17, the 
angular momentum measured by distant observers is en- 
hanced as compared to standard GR. We evaluate the 
correction to the angular momentum A J^s = Jcs ~ Jgr 
as a function of the coupling strength and find that it 
increases quadratically with the coupling parameter in 
the small-coupling regime, and asymptotes to a constant 
value A Jinax ~ {M / R)Jq^ in the large-coupling regime. 

These results are obtained using an analytic approxi- 
mation for constant-density nonrelativistic objects, and 
confirmed with a numerical solution for neutron stars de- 
scribed with realistic equations of state. The black hole 
solution is computed numerically and checked against 
known analytic solutions in the small coupling regime. 

Finally, using measurements of frame-dragging effects 
around the Earth, we set a weak but robust constraint 
to the CS characteristic lengthscale, S}/'^ < 10* km. Wc 
argue that this bound is the only current astrophysical 
constraint to the theory. 

This paper is organized as follows: in Sec. |ITj we re- 
view the theory of dynamical Chern-Simons gravity and 
define our notation. In Sec. Ill we lay out the general 



iting cases in Sec. IV as a well as an analytic solution 



formalism to compute the CS scalar field and the metric 
in CS gravity, for a slowly rotating object. We provide 
analytic expressions for the exterior solution in some lim- 



for the full spacetime for nonrelativistic constant-density 
stars. We present the results of our numerical computa- 
tions for realistic neutron stars in Sec. |V] and for black 
holes in Sec. 



VI 



We discuss constraints to the theory in 

SecEniand conclude in Sec. Ivml 

Throughout this paper we use geometric units G = c = 
1. We adopt the conventions of Ref. jTTj for the signature 
of the metric, Riemann and Einstein tensors. 



II. DYNAMICAL CHERN-SIMONS GRAVITY 

For a review on Chern-Simons modified gravity, we 
refer the reader to Ref. 2J. Here we simply recall the 
main equations and results and define our notation. 

We consider the following action defining the modified 
theory 




(1) 



In Eq. ([T]), the first term contains the standard Einstein- 
Hilbert action and the matter contribution with la- 
grangian density /Cmat- The second term is the CS mod- 
ification, which only depends on the dimensionfuU cou- 
pling constant £cs (which has dimensions of length; its 
relation to the parameter ^ of Ref. [7] is ^ = 167r^^g), and 
of course on the shape of the potential V for the dimen- 
sionless CS scalar field t9. The Pontryagin density RR 
is given by the contraction of the Riemann tensor Rap^v 
with its dual: 



RR 



~ 2 



(2) 



where e^'"^'^ is the four-dimensional Levi-Civita tensor 
(with convention e'^i^s = a right-handed orthonor- 

mal basis). The equations of motion resulting from the 
modified action are the following: 

(i) The modified Einstein field equations: 



(3) 



where G^i, is the Einstein tensor, T^'Jf* is the matter 
stress-energy tensor. 



^ The conversion from the notation used in Ref. 2^ (AY09) to 
that of the present work is given by i9 = \//3i?ayo9i IStt^^j, = 
= ^, = I3Vayo9{'S/VP)- The conversion from the 

notation of Ref. |4| (who define the Riemann tensor with an 
opposite sign) is given by =1/3. 
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is a four-dimensional generalization of the Cotton- York 
tensor, and 

= V^79V,i9 - ^.g,,.V„W"z? - g^uVi-d) (5) 
is the stress-energy tensor associated with the scalar field 

(a) The evolution equation for the scalar field "d: 

£^ 

ad ^ -^RR+V'i'd), (6) 

where □ = 5^''V^V,y is the usual covariant 
d'Alembertian operator. 

In general, the fundamental theory from which CS 
gravity arises should predict a shape for the potential 
This would introduce an additional free param- 
eter (or several parameters) in addition to i^^ . Since 
our goal here is mainly to study the effect of the cou- 
pling strength ^^g, we follow previous works and assume 
= for simplicity. 

III. SLOWLY-ROTATING RELATIVISTIC 
STARS IN CHERN-SIMONS GRAVITY 

A. Metric 

We consider a stationary, axially symmetric system. 
This means that there exist a time coordinate x*^ — t 
and an angular coordinate = (j) such that the metric 
components do not depend on t nor (j)] if x^ and x'^ are 
the two remaining spatial coordinates, we have 

g^Lu ^ 9,,u{x^,x^)- (7) 

We moreover assume that the line element ds^ — 
gfj^^dx^dx'^ satisfies an additional discrete symmetry: we 
suppose that the effect of inversion of the azimuthal angle 
is identical to that of inverting time, </> — > —0 t ^ —t, 
i.e. 

ds'^{t, x^jx'^, —(f); dt, dx^,dx^, —dcf) (8) 
— ds'^{~t,x^ , x'^ , <j)', —dt,dx^ ,dx'^ , d(j)). (9) 

Note that one could clearly not have made this assump- 
tion in nondynamical CS gravity where the externally 
prescribed scalar field defines the fiow of time (for the 
"canonical" choice of scalar field). The validity of this 
assumption in dynamical CS gravity will be justified a 
posteriori by checking that solutions do exist with this 
discrete symmetry (however, it does not guarantee the 
uniqueness of the solutions found). With this assump- 
tion, we find thai[^ 

3txi = gtx^ = 30x1 = g^^2 = 0. (10) 



^ In GR, one does not need the additional assumption |9| to obtain 
ITot , see proof in Section 7.1 of Ref. Il2j. The proof of Ref. 112] 



By an appropriate choice of coordinate transformations 
of the form i* = /'(x^,x^) for i = 1,2, one may rewrite 
the line element in the form [T3] 

-t- r^H{r, 6) [dO^ + (sin e fdif^) 

- (sin 6l)2w(r, 6')did0, (11) 

where we have removed the tildes for clarity and denoted 
x^ = r, = 0. 

B. Stress-energy tensor 

We consider an ideal fiuid in solid rotation with angular 
rate dip/dt = H., i.e. such that the 4- velocity is of the form 

= {u*, u'', w", w"^} = w°(r, 9){1, 0, 0, n}. (12) 

To evaluate u", we use the normalization condition 
g^.yuf'u'' = -1, i.e. 

= (e^* + 2r2(sin0)2cjl] + r'^ H{?,inefn'y^''' . (13) 

The stress-energy tensor of a perfect fluid is then ob- 
tained from 

T^''' = {p + p)u''u'' +pg>'\ (14) 

where p and p are the density and pressure of the fluid 
as measured in its local rest frame. 



C. Slow-rotation expansion 

We consider a staiEJof radius R in solid rotation with 
angular velocity f2. The star is in slow rotation if U,R <C 
1. In what follows we shall expand the metric functions 
$, A, H, bj, gtr and gtg as well as the density and pressure 
fields in powers of VlR. 

Under the change of coordinates (j)' = — (/), the angu- 
lar velocity changes sign, W = -Vl, gtt, grr, gee, g<i,<p and 
therefore <&, A, H remain unchanged, while gt^ and there- 
fore Lu change sign. Similarly, considerations of the trans- 
formation properties of the stress-energy tensor show 
that p and p remain unchanged. Since the field equa- 
tion expressed in the t, r, 6, (j)' coordinate system are the 
same as those in the r, 0, (j) coordinate system for a fluid 
with angular velocity —D, (because there is no functional 



makes use of the GR Einstein field equation, however, and can 
therefore not be carried over to CS gravity. One could also retain 
the a priori non-zero gtr and g^g metric components and show 
that the CS field equations lead them to vanish. We have not 
chosen this path in order ot avoid excessive algebra. 
In this paper we use the word 'star' to refer to a general (non- 
empty) astrophysical object. 
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dependence on (p), we conclude that <I>, A, H, p and p are 
even functions of whereas w is an odd function of f2. 

In the slow-rotation approximation, we only keep the 
lowest-order contribution to each function. From the 
previous discussion, this means that we keep only the 
0(1) terms in $,A, i/, p and p and the ©(fii?) term in 
u. The next-order contributions are of relative order 
fii^jM^) X (riR)^ in CS gravity, where / is an un- 
known function of the coupling strength. It is beyond 
the scope of this paper to evaluate /, but we emphasize 
that our slow-rotation approximation is only valid as long 
as / X (QR)^ < 1. 

To lowest order, the velocity of the fluid becomes 



e-^*{l,0,0,17}. 



(15) 



Moreover, the Pontryagin density is of order 0{ftR) (we 
give the exact expression below). Therefore, the scalar 
field is also of order 0{nR) [from Eq.^ with V = 
0]. As a consequence, the components of the scalar field 
stress-energy tensor T^f^ [Eq. ([s])] are of order 0{nR)'^ 
and need not be considered at this level of approximation, 
i.e. the energy of the scalar field does not curve space- 
time to first order in the rotation rate. 



D. Stellar structure equations for $, A and H 

The leading-order contribution to the metric compo- 
nents which are even in fl can be obtained in the non- 
rotating case, = 0. In that case the system is spher- 
ically symmetric and the Pontryagin density vanishes. 
The scalar field being unsourced, the Einstein field equa- 
tions are identical to those of standard GR. By spherical 
symmetry, $,A, 7?, p and p are functions of r only. By 
rescaling the radial coordinate, one can moreover set H 
to unity. The functions $(r) and A(r) then satisfy the 
usual relativistic stellar structure equations |11) : 



A(r) 



4ln(l 



with 



2m(r) 



m,{r) = / 'inr''^ p{r')d'"' 



(16) 



(17) 



being the mass enclosed within a radius r. Inside the 
star, the function <I>(r) is the solution of the equation 

m + Airr^p 



$'(r) 



(18) 



r{r — 2m) 

where the pressure p{r) must satisfy the equation of hy- 
drostatic equilibrium 

+ (19) 

r(r — Zm) 

Outside the star, we have p = p = 0, and 

1 . / 2M\ , , 

1 1 , (20) 



$(r) = -Mr) - In 

2 \ r J 

where AI is the total mass of the star. 



E. Gravitomagnetic sector 

1. General equations 

Let US now consider the t(f> component of the metric. 
Evaluating the t(j) component of the modified field equa- 
tion to lowest order in fii?, we find, in agreement with 
Ref. [To] and Refs. [71 (HI outside the star: 



2m\ d'^Lu 
r J dr"^ 



3cotan6' 



1 - 

m 



2m 
r 



tt{p 
167r(p + 



m 
r 



An 



-pr 



l-p)r2 

p){n 

dd 
do 



duj 
dr 



(21) 



In addition, the evolution equation for the CS scalar field 
becomes, in the time-independent case and to first order 
in nR: 



/ 27n\ d^-d 2 / ra ^ , ^ ^\ 



dr 



cote 



dd 



RR, 



(22) 



where, to lowest order, the Pontryagin density is given 
by 



RR 



-24- 



-(*+A) 



sm ( 



TO 



— p I — I sm"' 9 



89 



dr 



(23) 

Finally, we also note that all other non-diagonal com- 
ponents of the Einstein, "C" and stress-energy tensors 
vanish. Therefore the metric ( |11| is indeed a valid solu- 
tion of the modified Einstein field equations provided we 
solve for <&, A and w as described above. 



2. Multipole expansion 

Following Ref. [8], we decompose 'd{r,9) on the basis 
of Legendre polynomials and w(r, 9) on the basis of their 
derivatives: 



d{r,9) = Y,Mr)Pi{cost 

1=0 

00 

c^(r,^?) = ^L.,(r)P/(cos( 



(24) 
(25) 



1=1 



Eqs. (21 1 and (22) can then be rewritten as an infinite 



set of coupled differential equations for the coefhcients 
uji and I?;. Each pair of equations is independent in the 
sense that the l-th. multipole of uj only couples to the l-th 
multipole of i? [8]. The only sourced multipole is Z = 1, 
and therefore = and = for ^ ^ 1 (this follows 
from the system being second order in and "di , and the 
requirement that "d and uj vanish at large radii and be 
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bounded at the stellar center). The functions di{r) and 
tL>i{r) = a;(r) are the solutions of the following system 



1 

r J r 
+IQ'k{p + p){VL-uj) 



2m 



967r%e*-^ 



dr 



m Att 



(26) 



and 



/ 2m \ 2 / TO , s '?\ n, 2 „ 



12ee-(*+^) 



r \ 


TO 


r 


( m 


An 



p ^ 



(27) 



One can decompose uj into two pieces: a part that would 
be present in standard GR, cjg Ri that can be obtained 
by setting = in Eq. (26), and a correction (not 
necessarily small) Aw^s = lo — ujq^, which is the solution 
of Eq. ([26]) with rj = 0. Note that it is the full uj = 
WGR + Awcs that sources the CS scalar field in Eq. (27 1. 



Before proceeding further, we shall discuss boundary 
conditions at the surface of the star. First, the function 
Lj must be continuous at the boundary of the star. In 
addition, integrating Eq. ( 26 1 between R — e and R + e 



gives us a jump condition for the derivative of uj at the 
stellar surface: 



u'{R 



uj'{R- 



967r£2^ 47r 
~R2 §- 



piR~e)d,iR). (28) 



For neutron stars, the surface density is nearly vanishing 
[in the sense that p{R — e) ^ p(0)]i and uj' is (nearly) 
continuous at the surface. For constant density objects, 
however, there is a jump in the derivative of uj at the 
surface. Finally, inspection of Eq. (27) shows that all 



coefficients are bounded (although potentially discontin- 
uous), and therefore both "di and must be continuous 
at the stellar surface, in all cases. 



3. Simplification outside the star 



Outside the star, p — p ^ and Eq. ( 26 1 can be inte- 
grated once and simplified to 



6J 



(29) 



where J is a constant of integration. J is also the total 
angular momentum of the star as measured by observers 
in the asymptotically flat far zone, as we shall discuss in 
Sec. |IIIF| We write J = Jgr + A Jcs, where Jqr is the 
value of J in standard GR and AJ^s is the correction 
(not necessarily perturbative) that arises in GS gravity. 



Using Eq. (29) into Eq. (27), we obtain an equation for 
z9i alone: 

^ 2Af\ „ 2/ 

r J ^ r \ r J ^ yr^ J ^ 

^~12ll^, (30) 

where we have defined the dimensionless coupling 
strength C,, to be used repeatedly in the remainder of 
this paper: 



= 1287r 



(31) 



Physically, C}^'^ is of the order of the ratio of the CS 
lengthscale to the dynamical timescale of the system. It 
is important to notice that the definition of C, depends on 
the system considered through its average density. 
Equation ( 30 ) does not have an analytic solution in the 



general case, but it does have one in the two limiting cases 
of small coupling and nonrelativistic stars, discussed in 
Sec.llVl 



Equations (29) and (30 1 are also valid for a slowly- 
rotating black hole. Because a black hole does not have 
a surface of discontinuity but has a horizon, the boundary 
conditions for the scalar field and the metric are different 
than for a star. We shall discuss the black hole solution 
in Sec. IVD 



F. A note on the angular momentum and moment 
of inertia 

There are several possible definitions for angular mo- 
mentum in GR, and we therefore specify the definition 
that we use here. 

The angular momentum J of a body can be defined 
from the asymptotic behavior of the metric far outside 
the source [TTIITI] : 



As' 



1 



2M 

r 



o 



2M 

r 



O 



dtdx' 



5ij dx'^dx-' 



[l-%]dxMa;J'. 



(32) 



In spherical polar coordinates, and using the notation of 
Eq. (11), this corresponds to 



2J 



[l + 0(l/r)] 



(33) 



We therefore see that the constant of integration J in 
Eq. (29) corresponds to the total angular momentum of 



the star or black hole as measured by distant observers. 
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The moment of inertia is an ill-defined quantity for rel- 
ativistic systems, as in general the angular momentum of 
a body in solid rotation may not scale linearly with the 
angular velocity. However, in the slow-rotation approxi- 
mation, the angular momentum does scale linearly with 
fi, to first order. We can therefore define the relativis- 
tic generalization of the Newtonian moment of inertia by 



hm 



(34) 



In general the moment of inertia will depend on the mass 
and equation of state of the considered object. In Chern- 
Simons gravity, it will also depend on the coupling con- 
stant 

In GR, there exists a simple integral formula for the 
moment of inertia |15j : 



'GR 



Stt 



dr. 



(35) 



This formula is a priori valid only in GR, and does not 
necessarily hold in modified gravity theories. Equation 
( 35 1 is indeed a local definition, whereas we have defined 



angular momentum from its imprint on the spacetime in 
the far-field. 



It turns out, however, that Eq. (351 still holds in CS 
gravity, due to the particular form of Eq. (26 1. To see 



this, let us first rewrite Eq. (26 I in the form 



1 d 

r^ dr 

= 96n 



(r4e-(*+^'w') + 16^(p + p)e^-* {n - w) 



dr 



m 
r 



An 



-pr' z9i 



(36) 



Multiplying Eq. (36) by r , integrating form to i? — 
e, and using the jump condition for uj' at the surface 
of the star, Eq. (28), in conjunction with Eq. (29 1, we 
recover Eq. ( 35 ). We emphasize that this integral formula 



is not a definition of the moment of inertia and would not 
necessarily be valid in other modified gravity theories. 



IV. ANALYTIC APPROXIMATE SOLUTIONS 

Before tackling the full numerical solution of the prob- 
lem, we give a few analytic results in some simple cases. 



Analytic exterior solution in the small coupling 
limit 



If C ^ 1, the differential equation for -di outside the 



star Eq. ( 30 1 has an analytic solution [71 |B] , valid up to 



corrections of relative order ©(C^ 
i9i(C< 1,^ > R) 



5 J ^cs 

8 AP A 'P 



NP _^ 2M3 



18 Af* 



+C2(r--M) + 3Ci 



1 



\2M 



1 



ln(l - 2M/r 



where Ci and C2 are constants of integration. The re- 
quirement that § remains finite at infinity implies that 
C2 = 0. In Refs. [7j and where a rotating black hole 
was studied, Ci was also (rightly) set to zero so d (as well 
as w) remains finite at the horizon. In our case, however, 
since r > 2M outside the star, the homogeneous solution 
proportional to Ci is well behaved everywhere outside 
the star, and Ci 7^ a priori (in fact, we shall show 
that Ci sa 1 for nonrclativistic objects). The integration 
constant Ci must be determined from the continuity and 
smoothness requirements for ^ at the stellar boundary. 
This shows that the black hole solution cannot be used 
as the solution outside a star as was assumed in Ref. 1101. 



From Eq. (29), we then obtain uj outside the star (set- 
ting the additional integration constant to zero so that 
w(-Hoo) = 0): 



w(C< l,r > R) = 



2.J 



IOtt 



AP 



AP 12 



27 

loTs" 



15 M 
32T 



-Mn(l - 2M/r) 



15 AP 
3272" 

15 AP 

64 ^ 73" 



^AP 81 A/P 
?>~P ~ 16 7^ 

9AP 

47^ 



(38) 



It will be useful in what follows to write the asymptotic 
behavior of and a; at large radii r 3> M up to correc- 
tions of relative order Aijr (a fortiori, these results are 
also valid everywhere outside the source in the nonrcla- 
tivistic hmit i? > M): 



i9i(C< l,r> Af) 
a;(C< l,r>M) 



5 J ^cs 

8 AP AP 
2J 

jet 



(i-Ci) 



AP 



r 



IOtt 



AP 



(Ci -1 



aP 

^6 



32^ M5" 

5 r° 



64 M^l 



15 



Before going further, let us assess the differences of our 
results with those of Ref. [10], who also worked in the 
small coupling regime. First, as we pointed out previ- 
ously, Ref. Hn] set Ci = whereas we shall show below 
that Ci^l + 0{AI/R) + OiC). Therefore our solution 
for the scalar field is of order 0{M/R) smaller than that 
of Ref. [TD] at large distances. Second, and more impor- 
tantly, in Ref. [TD] , the parameter J used in Eq. ( 38 ) was 
set to Jgr- In reality, J = Jqr + ^^^CS is determined 
by imposing the continuity and jump conditions for lj at 
the surface of the star. In the small coupling regime. 



(41) 



,(39) 



(40) 



7 



and AJcs is in general nonzero. We therefore have 
Awes JC/f^ Eit large distances, instead of Awes 
{R^/A'P)JC/r^ in Ref. [TO]. Physically, this means that 
the CS correction translates into a change of angular mo- 
mentum as measured by distant observers, whereas the 
previous solution did not, strictly speaking, lead to any 
additional angular momentum. 



B. Analytic exterior solution in the nonrelativistic 

limit 



If M < R then the solution of Eq. ^ is 
J 



Mr >R:>M) 



'2t:CR^ 



1 ~ Do cosh 



-Di sinh 



CR^ 



(42) 



where Di and D2 are integration constants. Requiring 
d to be finite at large radii implies D2 = 1, whereas 
needs to be fixed using continuity conditions at the stellar 
boundary. 



From Eq. ( 29 1 , we then obtain uj outside the star 



(choosing the additional integration constant such that 
a;(+oo) = 0): 



w(r > i? > M) 



2J j fCR^' 

rh\ ""^ 



1 — cosh 



(43) 



If we Taylor-expand these solutions for ^ <C 1, we should 
recover the analytic solution obtained in Sec. IV A in 



the small-coupling limit, in the far-field limit, r M 
(we have just reversed the order in which the limits are 
taken). We obtain, up to corrections of relative order 



t?i(r->i?>M,C<l) 
oj{r > i?> M,C < 1) 



JDi 
\/2^r2 
2J 



(44) 



8V27r 



Jll,MD^ 



1287r J£1m2 



(45) 



Comparing to Eqs. (39) and (40), we see that the con- 



stants Ci and Di are related in the regime M R, ( ^ 
1: 



(46) 



Moreover, inspection of the second terms in Eqs. (39) and 
( 44 ) show that one must have 



Ci = l + 0{M/R)+OiC^). 



(47) 



C. Analytic solution in the nonrelativistic limit for 
a constant density object 

We have already obtained the general solution outside 
a nonrelativistic star in the previous section. We now 
need a solution inside the star to obtain the integration 
constants Di and J. 

Let us first consider Eq. (27) for di. For a constant 



density object, the right-hand-side vanishes, and in the 
non-relativistic limit, wc find 



Mr<R) = C3 



R 



C4 



r:' 



(48) 



Requiring "di to be finite at the center of the star, we 
set C4 = 0. We therefore have R^R - e) ^ di{R- e). 
The continuity and smoothness of "i^i at the boundary 
therefore imply that this relation is also satisfied at R+e. 
Imposing this condition with di given in Eq. (42 ) (where 
we recall that D2 = 1), we obtain, up to corrections of 
order AI/R, 



J 



C3 = , 

V27rCi?2 

Di — tanhC. 



1 



1 



cosh^ 



(49) 
(50) 



From the value of Di and Eq. (44) we infer the asymp- 
totic behavior of -di at large radii (r ^ R(^^/^): 



i?2 



^1 



J 



27ri?2 



tanh^. 



(51) 



We show the radial dependence of the scalar field di{r) 
for a nonrelativistic constant density object in Fig. [T] 
We see that in the linear regime (C < 1), the scalar field 
increases uniformly with C — for ^ ^ 1, i?i cx C^- In- 
creasing C further eventually leads to a damping of the 
scalar field near the surface of the star. In the far zone, 
we see from Eq. (51) that the asymptotic value of r'^'di 



plateaus to a constant value. 



Let us now consider Eq. (26 1 for to. Again, for a con 



stant density object, the right-hand-side vanishes. Using 
p = 3M/(47ri?'^), the equation satisfied by uj becomes, to 
first order in M/R: 

(52) 

which has the general solution, up to corrections of order 
{M/Rf: 



LJ = 9. + D-,\l 



6M 
5i?3' 



^ h _ iM,2 



(53) 



In order for w to be finite at the origin, we impose Z?4 — 0. 
Requiring uj to be continuous at the boundary and its 
derivative to satisfy the jump condition Eq. (28 1, wc find 
D3 = -| J/(Afi?2j and 



J = -MR^n 
5 



f - 



2M 



AM 

Jr 



tanh C 

c 



0{M/Rf 
(54) 
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FIG. 1. Chern-Simons scalar field i9i as a function of radius, 
for a constant-density nonrelativistic object. The curves are 
parametrized by the dimensionless coupling constant ^ de- 
fined in Eq. (31 1. For the ease of visualization, we have used 



dotted lines for C < 1 (the linear or quasi-linear regime, where 
i9i increases uniformly with (), and solid lines for C > 1 where 
non-linearity results into a damping of the scalar field in the 
near-zone and a plateauing of its asymptotic behavior at large 
radii. 



If we write J = Jgr + A Jcg (with the well-known result 
^GR = iMR'^n in the non-relativistic limit), we there- 
fore obtain, to lowest order in M /R, 

AJcs _ 4M tanhC' 
~ 5^ 



1 



Jgr 5i? V C 
In the small coupling limit (C ^ 1), we find 
AJcs _ 4AfC2 5127r£4^M3 



(55) 



(56) 



Jgr C«i 5i? 3 15 i?^ 

Interestingly, for large values of the coupling constant, 
the Chern-Simons correction saturates 
AJcs _ 4M 
Jgr C>i 5i? ■ 



(57) 



We show the function w(r) for a constant-density, non- 
relativistic object in Fig. [2j The effect of CS gravity 
is to (potentially strongly) decrease the gravitomagnetic 
potential near the surface, and enhance its asymptotic 
value in the far field r ^ RC^^"^ by a relative amount 
AJcs/Jgr < 4AJ/(5i?) « 1. 

In the next section we will see that the qualitative fea- 
tures of the analytic solution are recovered in the full 
numerical solution. 



NUMERICAL SOLUTION FOR REALISTIC 
NEUTRON STARS 



FIG. 2. Function ij{r) in CS gravity around a constant- 
density nonrelativistic object, for various values of the dimen- 
sionless coupling strength ( defined in Eq. (|31|), marked as la- 



bels. The dotted curve corresponds to standard GR (^ = 0). 
Because of the sharp boundary of a constant-density object, 
this function is continuous but not smooth, and the jump in 
its derivative at the surface is given by Eq. ( 28 1 . All curves 
asymptote to ui ~ 2J/r^ for r S> R(^^^^ but differ significantly 
close to the surface. The vertical dashed lines show the lo- 
cations probed by the LAGEOS [T^ and GPB ^ satellites, 
which orbit the Earth at an altitude of 6000 km and 640 km, 
respectively, and are used to set constraints to the CS length- 
scale in Sec. IVIII 



A. Equations of state 



Matter at nuclear densities has complex properties, 
and the EOS of neutron stars is not very well known. 
There exist two EOSs that are widely used in astrophys- 
ical simulations (see for example Ref. [T5] for a discus- 
sion): the Lattimer-Swesty (LS) EOS [T^ and the Shen 
et al EOS [101 HI]- We will use the LS EOS with nuclear 
incompressibility Kq = 220 MeV (hereafter, LS220) and 
the Shen et al EOS (hereafter Shen). We use the EOS 
routines of O'Connor and Ot^ [TH] to solve the relativis- 
tic stellar structure equations. We show the resulting 
mass-radius relations in Fig. |3] 

We also show the moment of inertia in standard GR 
as a function of neutron star mass and equation of state 
in Fig. |4] To eval uate the moment of inertia, we have 
integrated Eq. (36) in the absence of CS coupling (£^g = 



0) with a second-order implicit Euler method. 



In the previous section we have given analytic solu- 
tions for the coupled CS scalar field-metric system in the 
nonrelativistic limit, for constant density objects. In the 
present section we provide full numerical solutions for 
neutron stars described by realistic equations of state. 



* Available at |http; / / stellarcoUapse.org/microphysicsl we thank 
Evan O'Connor for providing tabulated solutions to the rela- 
tivistic equations of stellar structure. 
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FIG. 3. Neutron star mass-radius relation for the two equa- 
tions of state used in this work: Lattimer-Swesty (LS220) li9j 
and Shen et al PUIET] . 



are reversed, i.e. we have 



lim r^uj(r) = 2 J, lim —[r^uj{r)] ^ 0, (60) 

r— ^oo r— ^-oo qt* 

lim r^di{r) ^ d^, lim -^[r^Mr)] = 0, (61) 

r—>oo r— >-oo (XT 



where d4 and J are integration constants, and J — Jgr + 
A Jcs can physically be interpreted as the total perceived 
angular momentum of the system. 

In addition to these boundary conditions at r = and 
+00, the functions uj, uj' , "di and must all be continuous 
at the boundary of the star (in principle there is a jump in 
the derivative of ui, see Eq. (28 1, but for realistic neutron 



stars the surface density is 7 orders of magnitude than the 
mean density, so uj' is continuous up to small corrections) . 



0.0 L_ 

0.0 



LS220 
Shen 



2.0 



M/M., 



FIG. 4. Moment of inertia as a function of neutron star mass 
in standard GR, for the two equations of state used in this 
work: Lattimer-Swesty (LS220) il9i and Shen et al [^1^ . 



2. Shooting method 



We integrate the coupled system given by Eqs. (21) 



and (27) with a second-order implicit Euler method from 
r = to i?, and from r = +oo down tor = R (specifically, 
we consider the functions J)(u) = ^aj(l/u) and 'di{u) = 
-■di{l/u) from u = to 1/R; these functions vanish at 
u = and their derivatives at w = are proportional 
to the constants J and ^4). The linearity of the system 
allows to find the appropriate constants of integration 
at r = and -foo with a shooting method by requiring 
continuity and smoothness at r — R. 



Results 



B. Numerical solution of the coupled CS equations 

1. Boundary conditions 



Inspection of the system given by Eqs. (|21[) and (27) 



near the origin shows that there is one well behaved 
solution for w with w(r) ^ ci + 0{r^) and a diver- 

r— ^0 

gent (hence unphysical) solution with w(r) ^ C2/r'^. 

r— ^0 

Similarly, there is one well-behaved solution for with 
i?i(r) ^ c^r and a divergent, unphysical solution with 

r->0 



^1 



r->0 



c^/r . The physically allowed boundary condi- 



tions at the origin are therefore 

w(0)=ci, cj'(0)=0, 
^i(0)=0, M0) = C3, 



(58) 
(59) 



where ci and C3 are integration constants to be deter- 
mined. 

The mathematically allowed asymptotic behaviors at 
infinity are similar, but the physically relevant solutions 



In this section we illustrate our numerical results in 
several figures, and compare them to our analytic solu- 
tion of Sec. IrTcI 



1 . Scalar field 

First, in Fig. [s] we show the scalar field 'i?i(r) as a 
function of radius, for several values of the dimensionless 
CS coupling parameter (. A qualitative difference of this 
work with that of Ref. [lOj is that we have properly en- 
forced the continuity and the smoothness of the CS scalar 
field at the stellar boundary (see Fig. 4 of Ref. [TU] for 
a comparison). In the linear regime C < 1, t9i increases 
uniformly with For C ?i 1; more complex nonlinear be- 
haviors appear (the system solved being equivalent to a 
fourth-order ODE, such behaviors can be expected), but 
the overall effect is to suppress -d near the stellar surface. 
The far-field behavior of -d plateaus to an asymptotic 
limit, as we found in our analytic approximation. 
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FIG. 5. Radial dependence of the Chern-Simons scalar field 
I?! (r) for a 1 Mq neutron star described by the Shen et al EOS 
(the LS220 EOS shows a very similar behavior and we have 
not plotted it here for more clarity). The curves are lab elled 
by the dimensionless coupling strength C, defined in Eq. ( |31[ ). 
For C, the scalar field increases uniformly over the whole 
range r > 0. For C li the scalar field gets suppressed in the 
vicinity of the star; in the far field, lim r^i?i(r) asymptotes 

r — ^oc 

to a constant value for large values of ^. We have shown the 
two regimes C ^ 1 ^^^d C > 1 with different line styles for 
more visual clarity. 



2. Gravitomagnetic sector and moment of inertia 

In Fig. [6j we show the gravitomagnetic sector of the 
metric through the function uj{r). An essential differ- 
ence of our solution with previous work is that we have 
Awcs ~ 2AJcs/'"'^ at large radii, instead of Awcs oc 
l/r^. 

In Fig. [7] we show the relative change in moment of 
inertia induced by the CS modification, A/cs/-^gRi as a 
function of the dimensionless coupling strength ^. We re- 
call that the moment of inertia is defined as a / = J/^, 
where J — \ lim \r'^uj{r)\. We also plot the result of 

our analytic approximation for a nonrelativistic constant 
density star. Although the overall normalization is off by 
nearly an order of magnitude, we see that the trends pre- 
dicted by our simple approximation are indeed recovered 
in the full numerical result. In the limit of small CS cou- 
pling, we have A/cs/-^gr oc oc i'^^. For large values of 
the dimensionless coupling parameter, the CS correction 
to the moment of inertia asymptotes to a constant value. 

We show the dependence of A/cs /^gr on neutron star 
compactness in the small coupling limit in Fig. [Sj Again, 
our very simple analytic result is underestimating the ex- 
act numerical result, but the trend A/qs/^gr oc (M/R)^ 
is roughly respected for realistic neutron star masses. 

Finally, in Fig. [9] we show the extrapolated asymp- 
totic value of A/csZ-^gr in the limit of large coupling 
parameter, as a function of neutron star compactness. 
It depends approximately linearly on neutron star com- 
pactness, as we found in our analytic solution. To obtain 



FIG. 6. Radial dependence of the normalized metric coeffi- 
cient uj{r)/fl, for a 1 Mq neutron star described by the Shen 
et al EOS (the LS220 EOS shows a very similar behavior and 
the corresponding result was not plotted here). The curves 
are labelled by the dimensionless coupling strength defined 
in Eq. (31|. For C ^ 1 frame dragging becomes strongly sup- 
pressed in the vicinity of the star. At large radii uj ~ 2J /r^ , 
where J = Jqr -|- A Jcs is enhanced with respect to the value 
in GR. 



this value we have computed A/cs/^gr(C) for several 
values of the coupling strength up to ^ « 10 and fit- 
ted A/cs//gr = A/cs //gr loo (1 - ct/O — we could not 
go much beyond C « 10 as the homogeneous solutions 
contain terms of order exp(^) that quickly become very 
large and lead to large numerical errors, see equations of 
Sec UVBl 

Our numerical results are therefore in very good qual- 
itative agreement with our analytic approximation for 
constant-density nonrelativistic objects. Given that a 
neutron star density profile is far from constant and that 
neutron stars are compact, it is not surprising that the 
quantitative agreement is not perfect. We can under- 
stand why the analytic approximation systematically un- 
derestimates the correct result. A realistic neutron star 
is much more centrally concentrated than a constant den- 
sity object. The "effective" radius of the star (defined, 
for example, as the radius containing 75% of the mass), 
is always significantly smaller than the actual stellar ra- 
dius. A neutron star is therefore more compact in its 
central regions than if it had a constant density. For ex- 
ample, the mean density of a 1 Mq neutron star inside 
the sphere containing 75% of its mass is about 70% larger 
than its overall mean density. 



VI. SLOWLY-ROTATING BLACK HOLE IN 
CHERN-SIMONS GRAVITY 



In this section we compute the solution for the scalar 
field and the metric of a slowly-rotating black hole in CS 
gravity, for an arbitrary coupling constant. 

The equation satisfied by the CS scalar field for a 
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FIG. 7. Change in the moment of inertia A/cs/^gr induced 
by the CS modification as a function of the CS coupling 
strength for the two EOSs considered. We plot A/csZ-^gr 
for the values M/R = 0.1 and 0.15. The dashed lines show 
the analytic result derived in Sec. |IV C| for a constant den- 
sity nonrelativistic star (also for M/R = 0.1 and 0.15 from 
bottom to top). 
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FIG. 8. Relative change in the moment of inertia A7cs in- 
duced by the CS modification in the small coupling limit, 
as a function of neutron star compactness M/R, for the 
two EOSs considered. We show the dimensionless quantity 
limj_»o R"^ / its^Ics / Igk- The dashed line shows the analytic 
result derived in Sec. |IVCl for a constant density nonrelativis- 
tic star. The amplitude is off by approximately an order of 
magnitude but the overall behavior is relatively well repro- 
duced. 



slowly- rotating black hole is the same as that outside 
a slowly-rotating star, Eq. (30 1. Using x = r/(2M), and 



defining the dimensionless coupling parameter in analogy 
with Eq. (pTl) with the substitution R^2M 



a2 
5BH 



27r- 
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FIG. 9. Asymptotic value of the CS change in moment of 
inertia, max(AJcs/^GR), as a function of neutron star com- 
pactness. Also shown is our analytic approximation for non- 
relativistic constant density objects. 



the equation for the scalar field becomes 



1 



9Cbh 



9 



J 



^1 



(63) 



where dots denote differentiation with respect to x. 

There are two qualitative differences between the black 
hole case and the stellar case. 

Firstly, in the case of a slowly-rotating black hole, the 
angular momentum J is just a given parameter (as well as 
the black hole mass), independently of the gravity theory 
chosen. Of course, J is the angular momentum of the 
object that collapsed into a black hole, in which case its 
value does depend on the gravity theory chosen; once 
the black hole is formed, however, there is no way to 
disentangle Jgr from Jcs- 

Secondly, the boundary conditions are different than in 
the case of a star. For a black hole, the physical solution 
for d must be continuous and smooth at the horizon r — 
2M . Since there is a homogeneous solution that behaves 
as d ^ ln(l — 2M/r) near the horizon [see Eq. (37)], we 
must chose the boundary conditions in a way to avoid 
such a logarithmic divergence. 

To numerically solve Eq. (63), we start near the singu- 



lar point X — 1 with initial conditions 



^i(l 



C, 



i?i(l + e) = - 



9 ^csJ 

4 A/5 



(2 + gcin) c, 



(64) 



where C is a constant to be determined and the second 
equation follows from Eq. (63) evaluated at x = 1. We 



then integrate Eq. (63) outward to some radius xq > 1. 



We also integrate Eq. (63) from +oo to xq, with initial 
conditions lim [x'^-di{x)\= D, where D is a constant. 

We then adjust C and D such that i9i is continuous and 
smooth at the junction radius xq. Finally, w is obtained 



M4' 



(62) from Eq. ^ 
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FIG. 10. Chern-Simons scalar field outside the horizon of a 
slowly rotating black hole, for various values of the dimension- 
less coupling parameter C, = ^/2tv£'^^/M'^ . For better clarity, 
we have shown the small-coupling regime C ^ 1 with dotted 
lines and the nonlinear regime C > 1 with solid lines. 



VII. CONSTRAINTS TO THE THEORY 

To lowest order in the rotation rate, the Chern-Simons 
modification to general relativity only affects the gravit- 
omagnetic sector of the metric — for large enough values 
of the coupling strength, however, it is possible that the 
prefactor of terms of order (Ril)^ become significant in 
the gravitoelectric sector, but we have not addressed this 
issue in this paper. 

Tests of the theory will therefore rely on frame- 
dragging effects. The CS modification changes the t(j) 
metric component in two ways. First, it enhances the 
effective angular momentum J — Jgr + AJcs7 defined 
as J = 5 lim r^Lj. Secondly, it modifies the near-zone 

r— f CO 

(r ^ R) gravitomagnetic sector around the spinning ob- 
ject. We will use the latter effect to set a constraint on 
^cs with measurements of frame-dragging effects around 
the Earth, and the former effect to discuss potential 
constraints from the double-binary-pulsar system PSR 
J0737-3039 t22j. 



0.01 - 




FIG. 11. Metric coefficient a;(r) outside the horizon a slowly 
rotating black hole, for various values of the dimensionless CS 
coupling parameter = \/27r^^s/Af^. 



The solution scales linearly with the dimensionless pa- 
rameter J/M^ and otherwise depends on the single pa- 
rameter Cbh ■ We checked that in the small-coupling limit 
Cbh ^ 1 our numerical solution agrees with the analytic 
result of Refs. [71 [8]. We show the scalar field for sev- 
eral values of Cbh in Fig- [lO) and the metric coefficient 
uj in Fig. [TT] We see that the black hole solution ex- 
hibits the same features as the solution for a star: the 
scalar field first increases uniformly with Cbh and then 
starts being damped near the horizon in the nonlinear 
regime Cbh ^ 1- Frame-dragging effects become strongly 
screened next to the black hole horizon for large CS cou- 
pling strengths. More specifically, for Cbh ^ Ij screening 
occurs for 2M < r < 2MC 



Measurements of frame-dragging around the 
Earth 



The Earth is clearly a nonrelativistic object, with 
M<^/R^ « 7 X 10"^°. Moreover, the density of the Earth 
varies by no more than a factor of a few from center to 
edge, and the results of the analytic approximation pre- 
sented in Sec. |IV C| should provide fair estimates. The 
change in angular momentum A Jcs for an object so lit- 
tle compact as the Earth is completely negligible, with 
AJcs/'/gr < ^^/R < 10~^. However, the gravitomag- 
netic field close to the Earth's surface may be significantly 
modified for a large enough CS coupling constant, as one 
can see in Fig. [2j 

Gravity Probe B (GPB) has measured the gyroscopic 
precession due to frame-dragging to be in agreement with 
GR to an accuracy of 20% [T7] . The GPB satelhte is in 
orbit at an altitude of 640 km, which corresponds to a 
distance from the Earth's center r k, l.li?0. This means 
that we hav^E] 



c^(l.li?e) 



2J/(l.li?e) 



< 20%. 



(65) 



By requiring that the constraint Eq. ( 65 ) be satisfied and 
using Eq. ( [43| for w (with Di = tanhCe), we obtain Ce ^ 
1 (we emphasize that C being a dimensionless constant, it 



1/3 
BH • 



^ Frame-dragging is due to the gravitomagnetic field B = V X A, 
where the gravitomagnetic potential is A = Ir sin 6'aj(r)e^ in 
spherical polar coordinates. In principle the precession of gyro- 
scopes or the Lense-Thirring drag depend on oj and its derivative. 
For an approximate constraint as we give here, the detailed ex- 
pression is not crucial, though, the idea being that for C 5^ !> the 
gravitomagnetic field becomes significantly affected. 
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depends on the mass and radius of the system considered, 
hence the subscript ©). Translating this into a constraint 
to 4s, or, to use the notation of Ref. [7], the parameter 
^ = 167r£^g, we obtain: 



jl/4 ^ 



(167r)i/44s < 10^ km. 



(66) 



The LAGEOS and LAGEOS 2 sateUites have allowed 
to detect the Lense-Thirring effect for bodies orbiting 
the Earth with a 10% precision [1^]. This precision test 
of GR has allowed to set the first constraint on nondy- 
namical Chern-Simons gravity [4] (since then improved 
with measurements of the precession rate in the double- 
binary-pulsar system The LAGEOS satellites are 
at a distance r « 12000 km « 2Rm from the Earth cen- 



ter. Using again Eq. ( 65 1 with the appropriate distance 



and uncertainty, we obtain the constraint Ce ^ 2, which 
leads, up to a factor of ~ -\/2, to the same bound as 



Eq. (66) 



Even though these are relatively weak constraints, they 
have the advantage of being independent and very robust. 
Even with an error of 50% on either measurement, the 
constraint would be degraded by less than a factor of 3. 
Moreover, we will explain below that these are actually 
the only current astrophysical constraints (to our knowl- 
edge) of dynamical CS gravity. 



B. On future constraints from the 
double-binary-pulsar 

We have seen in the previous sections that CS gravity 
can only decrease frame-dragging effects near the surface 
of a spinning star, and only marginally enhance them 
in the far-field, even for an arbitrarily large coupling 
strength. To be able to test CS gravity, observational 
accuracy must therefore reach the level at which frame 
dragging effects can be measured (so one can estimate 
deviations from GR predictions). In other words, up- 
per bounds on frame dragging effects cannot be used to 
set any constraint on the CS coupling strength (unless 
the bounds are in fact lower than the GR predictions). 
Ref. [Tj derived the constraint ^^/^ < 1.5 x 10^ km from 
the measurement error on the periastron precession rate 
in the double-binary-pulsar system, which is much larger 
(by two orders of magnitude) than the GR-predicted pe- 
riastron precession rate due to spin-orbit coupling. This 
bound was derived by extrapolating the small-coupling 
result to large coupling strengths. In regard of the above 
discussion, we see that current measurements of the peri- 
astron precession rate in the double-binary pulsar in fact 
cannot be used to set any constraint on CS gravity. 

It has been suggested that long term measurements of 
the binary system PSR J0737-3039 [25] would lead to a 
10% determination of the moment of inertia of pulsar A 
[23l [24] , although this would likely take another 20 years 
of measurements |25j . Ref. have used this prediction 
to derive a potential future constraint on the CS cou- 
pling strength. From Fig. [vj we see that a 10% accuracy 



measurement of the moment of inertia would translate 
in a bound Cns ^ 1, which would lead to the constrainl]^ 
f < 25 km. Note that here we have assumed that pul- 
sar B lies in the far field of pulsar A, i.e. that the semi- 
major axis a « 4 X 10^ km is such that (a/R)^ ^ Cns- 
Using the constraint from Earth measurements Eq. ( 66 1 , 
we find that Cns is actually at most ~ 10^"^ ~ (a/R)^ and 
assuming the near-field limit would therefore not improve 
on the terrestrial constraint anyway. 

Whereas the prospect of such a strong constraint is 
enticing, we should keep in mind, first, that this is a 
very difficult measurement [25] . Furthermore, a second 
challenge is the inherent uncertainty in the GR-predicted 
moment of inertia, which varies by ^ 20 — 30% depending 
on the equation of state used (using the LS220 or Shen et 
al. EOSs). In fact, Ref. [H [25 initially suggested that 
the measurement of the moment of inertia could help nar- 
row down the space of allowed equations of state, since 
the variations in moment of inertia are larger than those 
in radius (by a factor of ~ 2 since / oc MR?), which is 
hard to measure in any case. We therefore conclude that 
constraining CS gravity from measurements of the mo- 
ment of inertia is likely to be very challenging, if possible 
at all. 



As a conclusion, the bound ( 66 ) that we have derived 



from measurements of frame-dragging around the Earth 
is, to our knowledge, the only current astrophysical con- 
straint to dynamical CS gravity. 



VIII. CONCLUSIONS 

We have presented the solution to the coupled system 
of field equations and the scalar field Poisson equation in 
dynamical CS gravity, for slowly rotating stars and black 
holes. We have shown that the black hole solution does 
not describe the spacetime outside a rotating star as was 
assumed in previous works. For the first time, we have 
provided a solution valid for an arbitrary CS coupling 
strength and not limited to the small-coupling regime. 
We have provided simple analytic solutions for nonrel- 
ativistic constant-density objects, as well as numerical 
solutions for realistic neutron stars and slowly-rotating 
black holes. 

Our solution shows two key features. First, frame- 
dragging effects are reduced with respect to standard GR 
near the surface of a rotating star or black hole. We 
have used this effect to set a robust constraint on the 
CS lengthscale from measurements of the Lense-Thirring 
drag and gyroscopic precession around the Earth, C^/'' < 



This constraint is of the same order as the one derived in Ref. 'lOJ , 
even though their metric solution actually does not contain any 
additional angular momentum, since they have &u) ~ at 
large radii. However, because they used the integral formulation 
( |35| for the moment of inertia, they did obtain the correct order 
of magnitude for the constraint. 
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10* km. Probes of the space-time close to the horizon of 
spinning black holes with orbits of objects passing nearby 
could potentially help constrain the theory further [5]. 
Secondly, the angular momentum of a rotating object (as 
perceived by observers in the far field) is enhanced. How- 
ever, this enhancement is at most AJcs/^gr ^ M/R 
and is therefore difficult to detect. 

In closing, let us mention some limitations of our work. 
Firstly, we recall that this paper only considers the slow- 
rotation approximation. Arbitrarily fast rotating stars or 
black hole solutions are significantly more complicated, 
even in GR, and we have not tackled this difficult problem 
here. Secondly, we have only studied the Chern-Simons 
extension to GR and have not considered any other ex- 
tension of quadratic (or higher) order in the Riemann 
tensor. Recently, Ref. [M| have considered compact stars 
in Einstein-Dilaton-Gauss-Bonnet gravity, which comple- 
ments the present work for quadratic extensions to GR. 
Finally, we have only studied stationary spacetimes, and 
have not performed any stability analysis which would 
require a full time-dependent study. It is therefore not 
clear that the solutions we have found are stable, and, in 
fact, it is likely that instabilities occur for large coupling 
strengths as the system solved is effectively fourth-order. 
It is also possible that higher-order terms in the "exact" 



theory of which the Chern-Simons extension is a trunca- 
tion render the solution stable. These questions would 
require a significantly more thorough analysis, which we 
defer to future work. 
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